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Abstract. The direct detection of dark matter is central to particle physics 
and cosmology. Current fashionable supersymmetric models provide a nat- 
ural dark matter candidate which is the lightest supersymmetric particle 
(LSP). Such models combined with fairly well understood physics like the 
quark substructure of the nucleon and the structure of the nucleus (form 
factor and/or spin response function), permit the evaluation of the event 
rate for LSP-nucleus elastic scattering. The thus obtained event rates are, 
however, very low or even undetectable. So it is imperative to exploit char- 
acteristic signatures, like the modulation effect, i.e. the dependence of the 
event rate on the earth's annual motion, and the directional rate, i.e its 
dependence on the direction of the recoiling nucleus. In this paper we do 
this using various velocity distributions, isothermal (symmetric as well as 
only axially asymmetric) and non isothermal (e.g. due to caustic rings). 



1. Introduction 

The consideration of exotic dark matter has become necessary in order to 
close the Universe [1]' [2]. Furthermore in in order to understand the large 
scale structure of the universe one has to consider matter made up of par- 
ticles which were non-relativistic at the time of freeze out, i.e cold dark 
matter (CDM). The COBE data [3] suggest that CDM is at least 60% 
[4]. Recent data from the High-z Supernova Search Team [5] and Super- 
nova Cosmology Project [6] ■ [7] suggest the necessity of the cosmological 
constant A. In fact the situation can adequately be described by a barionic 
component &b = 0.1 along with the exotic components Qcdm = 0.3 and 
= 0.6. In another analysis Turner [8] gives Q m = &cdm + = 0.4. 
Since the non exotic component cannot exceed 40% of the CDM [2]' [9], 
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there is room for the exotic WIMP's (Weakly Interacting Massive Parti- 
cles). In fact the DAMA experiment [10] has claimed the observation of 
one signal in direct detection of a WIMP, which with better statistics has 
subsequently been interpreted as a modulation signal [11]. 

The above developments are in line with particle physics considerations. 
Thus, in the currently favored supersymmetric (SUSY) extensions of the 
standard model, the most natural WIMP candidate is the LSP, i.e. the 
lightest supersymmetric particle. In the most favored scenarios the LSP 
can be simply described as a Majorana fermion, a linear combination of 
the neutral components of the gauginos and Higgsinos [2]' [12]~[23]. 

Since this particle is expected to be very massive, m x > 30GeU, and 
extremely non relativistic with average kinetic energy T < WOKeV, it can 
be directly detected [12]~[13] mainly via the recoiling of a nucleus (A,Z) 
in the elastic scattering process: 

X + (A Z) -> X + (A,Z)* (1) 

(x denotes the LSP). In order to compute the event rate one needs the 
following ingredients: 

1) An effective Lagrangian at the elementary particle (quark) level ob- 
tained in the framework of supersymmetry as described in Refs. [2], Bottino 
et al. [20], [23], Kane et al, Castano et al. and Arnowitt et al. Our own 
SUSY input parameters will appear elsewhere [19] : 2) A procedure in going 
from the quark to the nucleon level, i.e. a quark model for the nucleon. The 
results, for the scalar and the isoscalar axial couplings, depend crucially on 
the content of the nucleon in quarks other than u and d [14, 24]. 

3) Compute the relevant nuclear matrix elements [25]~[29]. 

Since the obtained rates are very low, one would like to be able to exploit 
the modulation of the event rates due to the earth's revolution around the 
sun. One also would like to know the directional rates, by observing the 
nucleus in a certain direction, which correlate with the motion of the sun 
around the center of the galaxy. 

The purpose of our present review is to focus on this last point along 
the lines suggested by our recent work [17, 18]. 

2. Expressions for the Unconvoluted Event Rates. 

Combining for results of the previous section we can write 

Leff = -^{(X!7 A 75Xi)A + (xiXi)J} (2) 

where 

Jx = N lx {f v + /i T 3 + /0 75 + /i 75 r 3 )iV , J = N{f s + fln)N (3) 



3 



We have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the Majorana nature of the LSP, Xi7 A Xi = (identically). 

With the above ingredients the differential cross section can be cast in 
the form 

Mu ' v) = 2(/I.)' [(gs + s " i ] p2{u) + s *" i " Fii( " )l (4) 



J spin 



spin > 



C*. = ^^[(^„(0)f §^ + 4^(0)^,(0)^1 + 0,(0))^] 

(7) 

_ (^/mjy) 2 2 fy A - 2Z 2 vp 2 1 2r? + l ( 2^ ) 

^ (1 + 4)2 1 f v A c ^ (2^6)Ml + r/) 2 (^ 2 ) J 

(9) 

<j l o = proton cross-section, i = S", spin, V given by: 

x° = a ° (fs") 2 ( sca l ar ) ) (the isovector scalar is negligible, i.e. a'p = a'^) 
af™ = * 3 U% + /A) 2 (spin) , a^ = a (/° + f v ) 2 (vector) 
where m p is the proton mass, r\ = m x /mj^A, and is the reduced mass 
and 

a = — (G F m N ) 2 ~ 0.77 x 10" 38 cm 2 (10) 

27T 

n = (7 2 6 2 /2 or Q = Q u, Q = A J ^ (11) 

where b is (the harmonic oscillator) size parameter, q is the momentum 

transfer to the nucleus, and Q is the energy transfer to the nucleus 

In the above expressions F(u) is the nuclear form factor and F pp >(u) arc 

the spin form factors [13] (p, p are isospin indices) 

Both form factors are normalized to one at u = 0. 

The nom-directional event rate is given by: 

dN p(0) m 

R = Rnon-dir = -J- = 1 a(u,v)\v\ (12) 

at m x AniN 
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Where p(0) = 0.3Gey/cm 3 is the LSP density in our vicinity and m is the 
detector mass The differential non-directional rate can be written as 

dR= P(0) A m da(u,v)\v\ (13) 
m x Am N 

where da(u, v) was given above. 

The directional differential rate [30] in the direction e is given by : 

dR dir = P(0) m v.eH(v.e) da(u, v) (14) 
m x AniN 27r 

where H is the Heaviside step function. The factor of l/2n is introduced, 
since we have chosen to normalize our results to the usual differential rate. 



3. Convolution of the Event Rate 

We have seen that the event rate for LSP-nucleus scattering depends on 
the relative LSP-target velocity. In this section we will examine the conse- 
quences of the earth's revolution around the sun (the effect of its rotation 
around its axis is expected to be negligible) i.e. the modulation effect. In 
practice this has been accomplished by assuming a consistent LSP velocity 
dispersion, such as a Maxwell distribution [2] or its extensions with only 
axial symmetry [31], already been discussed in the lierature [17, 18]. More 
recently other non-isothermal approaches, in the context velocity peaks and 
caustic rings, have been proposed, see e.g Sikivie et al [32]. 

In the last case one considers the late in-fall of dark matter into our 
galaxy producing flows of caustic rings. In particular the predictions of a 
self-similar model have been put forward as a possible scenario for dark 
matter density- velocity distribution, see e.g. Sikivie et al [32]. The implica- 
tions of such theoretical schemes and, in particular, the modulation effect 
are the subject of this section. 

Following Sikivie we will consider 2xiV caustic rings, (i,n) , i=(+.-) and 
n=l,2,...N (N=20 in the model of Sikivie et al), each of which contributes to 
the local density a fraction p n of the of the summed density p of each of the 
i = +, — . It contains WIMP like particles with velocity y n = (y nx , y ny , y nz ) 
in units of essentially the sun's velocity (vq = 220 Km/s), with respect 
to the galactic center. The z-axis is chosen in the direction of the disc's 
rotation, i.e. in the direction of the motion of the the sun, the y-axis is 
perpendicular to the plane of the galaxy and the x-axis is in the radial 
direction. We caution the reader that these axes are traditionally indicated 
by astronomers as e>0,e r ,e 2 ) respectively. The needed quantities are taken 
from the work of Sikivie et al [32] , via the definitions 
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Vn = v n/vO,Vnz = *W^0 = Vnz^'nx = v nr/v = VnxiVny = v nz/v = 

Vny, Pn = d n /p, p = £n=i d n and y n = {(ynz - l) 2 + vl y + yl x ] l/2 (for each 
flow +.-). This leads to a velocity distribution of the form: 

N 

f(v>) = ]T S(v - v y' n ) (15) 

n=l 

Since the axis of the ecliptic [13]. lies very close to the y, z plane the velocity 
of the earth around the sun is given by 

ve = Vo + vi = vo + vi(sinax — cosacos'yy + cosasin^/z) (16) 

where a is the phase of the earth's orbital motion, a = around second of 
June. 

One can now express the above distribution in the laboratory frame [18] 
by writing v = v + ve 

4. Expressions for the Differential Event Rate 

The mean value of the non-directional event rate of Eq. (13), is given by 

/dR\_p(0) m r da(u,v) ^ 

\ du I m x Am jy J ' du 

The above expression can be more conveniently written as 

(?) = ^v^(?) ( 18 ) 

\dul m x Am N v du 

where 

,dE, f \v\ „. .da(u,v) ,q , . 

l ^ = J 7m ,{v - VE, ^" v (19) 

There are now experiments under way aiming at measuring directional 
rates , i.e. the case in which the nucleus is observed in a certain direction. 
The rate will depend on the direction of observation, showing a strong 
correlation with the direction of the sun's motion. In a favorable situation 
the rate will merely be suppressed by a factor of 2tt relative to the non- 
directional rate. This is due to the fact that one does not now integrate over 
the azimuthal angle of the nuclear recoiling momentum. The directional rate 
will also show modulation due to the Earth's motion. 

The mean value of the directional differential event rate of Eq. (14), is 
defined by 



/dR\ p(0) m 1 r „ da(u,v) 

\ du I dir m x Am jy 2ir J ' E du 
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It can be more conveniently expressed as 

/dR\ _ p(0) m r-jr ■ dE 

\ du I dir m x AlTlN^ du dir 

where 

,dE. 1 f v.eH(v.e) „, .da(u,v) , ? . . 

< *: ) * = ^/^p^ /(t '' t ' £) -^ <^l • (22) 

We are not going to discuss the differential rates further. We will limit our 
discussion to the case of the total rates. 



5. The Total non-directional Event Rates 

Integrating the differential rate from u m in to u max we obtain for the total 
non-directional ratei. In the case of caustic rings we find: 

R = Rt^-[l-h(a,Q min )cosa)} (23) 

where 

R = ^-P-y/WPs + X spin + ^£y] (24) 

with Sj,i = S,V,spin have been defined above, see Eqs (5) - (8). Further- 
more a = [V2[i r bvo] _1 



Qmin ■ / Vesc t Vn \ • -, -. , - 

Umin = —^—,u max = mm(^-,maa:(^) , n = 1, 2, N) (25) 



Here y esc = with v escape = 625Km/s is the escape velocity from 

the galaxy. The modulation is described by the parameter h. Similarly 
integrating Eq. (19) we obtain for the total non-directional rate in our 
isothermal model as follows: 



R = Rt[{l + h(a,Q min )cosa)} (26) 

Now u max = v\ S cl ^ Note the difference of sign in the definition of the 
modulation amplitude h compared to Eq. (23). The modulation can be 
described in terms of the parameter h. The effect of folding with LSP 
velocity on the total rate is taken into account via the quantity t. 

The meaning of t is clear from the above discussion. We only like to 
stress that it is a common practice to extract the LSP nucleon cross section 
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from the the expected experimental event rates in order to compare with 
the SUSY predictions as a function of the LSP mass. 

6. The Total Directional Event Rates 

We will again examine the case of caustic rings and the isothermal models 
considered above. 

6.1. THE TOTAL DIRECTIONAL EVENT RATES IN NON-ISOTHERMAL 
MODELS 

We need distinguish distinguish the following cases: a) e has a component 
in the sun's direction of motion, i.e. < 9 < 7r/2, labeled by i=u (up), b) 
Detection in the direction specified by tt/2 < 9 < ir, labeled by i=d (down). 



Thus : 



1. In the first quadrant (azimuthal angle < </> < 7r/2). 




(27) 




+ {r l y + cosah l 2 (u) + -^-(\cosa\ — cosa))\e y .e 
U 

+ (r* + sinah\ + -^-{\sina\ + sina))\e x .e\} 
3. In the third quadrant (azimuthal angle it < (p < 3ir/2). 



(28) 



Ri dir = R^^-K r i- c osa h\)e z .e 




+ {r l y — cosah^u) H — ^ — (|cosa| — cosa))\e y .e 

+ (r* + sinaHl + -^-(\sina\ + sina))\e x .e\] 
4. In the fourth quadrant (azimuthal angle 3ir/2 < (p < 2ir) 




(29) 
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+ (r* — cosah\ + -^(|cosa| — cosa)) le^.ej 
+ i r x ~ sinah\ + -^-(\sina\ — sina))\e x .e\] 



(30) 



6.2. THE TOTAL DIRECTIONAL EVENT RATES IN ISOTHERMAL 
MODELS 

In this case we find it convenient to consider the absolute value of the 
difference of the rates in the two opposite directions. Integrating Eq. (20) 
and restricting ourselves close to the axes we obtain 



Note that now the rate is normalized to t°/2 and not to t and that the 
modulation can be described in terms of three parameters hi, 1=1,2,3. 

7. Results and Discussion 

We will discuss the effects of folding with the LSP velocity combined with 
the nuclear form factor and specialized in the case of the nucleus 127 1, which 
is one of the most popular targets [10]' [33] [34] _ [3 7] for the enegy cutoffs: 
Q min = 0, 10, 20 KeV. 

In the case of the non isothermal model of Sikivie et al [32], the total 
rates are described in terms of the quantities t,r l x ,r y ,r l z for the unmodu- 
lated amplitude and h, h\, h\, h\, h l c , h\ i = u,d for the modulated one. (see 
Table 1). Out of this list list only t and h enter the non-directional rate. We 
notice that the usual modulation amplitude h is smaller than the one aris- 
ing in isothermal models [17, 18]. As expected, the parameter t decreases 
as the reduced mass increases. 

In the case of isothermal models we will limit ourselves to the discus- 
sion of the directional rates. If the direction of observation is close to the 
coordinate axes, the rate is described in terms of the three quantities to 
and hi, i = 1,2,3 (see Eq. (31)). These are shown in tables 2-4 (for the 
differential rate see our previous work [17, 18]. 

8. Conclusions 

We have studied the parameters describing the event rates for direct detec- 
tion of SUSY dark matter and, in particular, two characteristic signatures, 
which will aid the experimentalists in reducing background: The non direc- 
tional modulated event rates, which are correlated with the motion of the 



Rdir = -R (t°/47r) |(1 + hi(a,Q min )cosa)e z .e 

- h 2 (a, Qmin ) cosae y .e + h 3 (a, Q min ) sinae x .e 



(31) 
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TABLE 1. The quantities t and h entering the total non-directional rate in the 
case of the target 53 J 127 for various LSP masses and Q m in in KeV. Also shown 
are the quantities r l j,hj i — u,d and j = x,y,z,c,s, entering the directional 
rate for no energy cutoff. For definitions see text. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t 


0.0 


1.451 


1.072 


0.751 


0.477 


0.379 


0.303 


0.173 


h 


0.0 


0.022 


0.023 


0.024 


0.025 


0.026 


0.026 


0.026 


r u 

1 z 


0.0 


0.726 


0.737 


0.747 


0.757 


0.760 


0.761 


0.761 


' y 


0.0 


0.246 


0.231 


0.219 


0.211 


0.209 


0.208 


0.208 


rl 


0.0 


0.335 


0.351 


0.366 


0.377 


0.380 


0.381 


0.381 


K 


0.0 


0.026 


0.027 


0.028 


0.029 


0.029 


0.030 


0.030 


K 


0.0 


0.021 


0.021 


0.020 


0.020 


0.019 


0.019 


0.019 


K 


0.0 


0.041 


0.044 


0.046 


0.048 


0.048 


0.049 


0.049 


n c 


o o 


o n^fi 

U.Uou 


o o^s 


U.U4U 


U.U41 


n D/19 

U.U4Z 




O 0/10 


K 


0.0 


0.036 


0.024 


0.024 


0.023 


0.023 


0.022 


0.022 


4 


0.0 


0.274 


0.263 


0.253 


0.243 


0.240 


0.239 


0.239 


r d 
1 y 


0.0 


0.019 


0.011 


0.008 


0.007 


0.007 


0.007 


0.007 


r d 


0.0 


0.245 


0.243 


0.236 


0.227 


0.225 


0.223 


0.223 


h d z 


0.0 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


K 


0.0 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


h d x 


0.0 


0.022 


0.021 


0.021 


0.020 


0.020 


0.020 


0.020 


h d c 


0.0 


0.019 


0.018 


0.018 


0.017 


0.017 


0.017 


0.017 


h d s 


0.0 


0.001 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


t 


10.0 


0.000 


0.226 


0.356 


0.265 


0.224 


0.172 


0.098 


h 


10.0 


0.000 


0.013 


0.023 


0.025 


0.025 


0.026 


0.026 


t 


20.0 


0.000 


0.013 


0.126 


0.139 


0.116 


0.095 


0.054 


h 


20.0 


0.000 


0.005 


0.017 


0.024 


0.025 


0.026 


0.026 



Earth and the directional event rates, which are correlated with both the 
velocity of the sun and the velocity of the Earth. 

In the case of non isothermal models [32] from Table 1 we see that 
the maximun no longer occurs around June 2nd, but about six months 
later. The difference between the maximum and the minimum is about 4%, 
smaller than that predicted by the symmetric isothermal models [17, 18]. 
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TABLE 2. The quantities t°,hi and h m for A = in the case of the target 
53 1 127 for various LSP masses and Qmin in KeV (for definitions see text). 
Only the scalar contribution is considered. Note that in this case and /13 
are constants equal to 0.117 and 0.135 respectively. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t° 


0.0 


1.960 


1.355 


0.886 


0.552 


0.442 


0.360 


0.212 


hi 


0.0 


0.059 


0.048 


0.037 


0.029 


0.027 


0.025 


0.023 


t° 


10. 


0.000 


0.365 


0.383 


0.280 


0.233 


0.194 


0.119 


hi 


10. 


0.000 


0.086 


0.054 


0.038 


0.033 


0.030 


0.025 


t° 


20. 


0.000 


0.080 


0.153 


0.136 


0.11 


0.102 


0.065 


hi 


20. 


0.000 


0.123 


0.073 


0.048 


0.041 


0.036 


0.028 



In the case of the directional rate we found that the rates depend on the 
direction of observation. The biggest rates are obtained, If the observation 
is made close to the direction of the sun's motion. The directional rates 
are suppressed compared to the usual non-directional rates by the factor 
fred = k/(2tt). We find that k = r" ~ 0.7, if the observation is made in 
the sun's direction of motion, while k ~ 0.3 in the opposite direction. The 
modulation is a bit larger than in the non-directional case, but the largest 
value, 8%, is not obtainrd along the sun's direction of motion, but in the 
x-direction (galactocentric direction). 

In the case of the isothermal models we restricted our discussion to the 
directional event rates. The reduction factor is now given by the parameter 
fred = W(47r t) = k/(2tt). We find that k is around 0.6 for no asymmetry 
and around 0.7 for maximum asymmetry (A = 1.0). In other words it is not 
very different from the naively expected f re d = l/(2ir). The modulation 
of the directional rate increases with the asymmetry parameter A and it 
depends, of course, on the direction of observation. For Q m in = it can 
reach values up to 23%. Larger values, up to 35%, are possible for large 
values of Q m i n , but they occur at the expense of the total number of counts. 

In all cases our results indicate that t for large reduced mass deviates 
from unity. Thus, if one attempts to extract the LSP-nucleon cross section 
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TABLE 3. The same as in the previous table, but for the value of the asym- 
metry parameter A = 0.5. 











LSP 


mass 


in GeV 






Quantity 




10 


30 


50 


80 


100 


125 


250 


t 


0.0 


2.309 


1.682 


1.153 


0.737 


0.595 


0.485 


0.288 


hi 


0.0 


0.138 


0.128 


0.117 


0.108 


0.105 


0.103 


0.100 


h 2 


0.0 


0.139 


0.137 


0.135 


0.133 


0.133 


0.133 


0.132 


h 3 


0.0 


0.175 


0.171 


0.167 


0.165 


0.163 


0.162 


0.162 


t° 


10. 


0.000 


0.376 


0.468 


0.365 


0.308 


0.259 


0.160 


hi 


10. 


0.000 


0.174 


0.139 


0.120 


0.114 


0.110 


0.103 


h 2 


10. 


0.000 


0.145 


0.138 


0.135 


0.134 


0.134 


0.133 


h 3 


10. 


0.000 


0.188 


0.174 


0.167 


0.165 


0.164 


0.162 


t° 


20. 


0.000 


0.063 


0.170 


0.171 


0.153 


0.134 


0.087 


hi 


20. 


0.000 


0.216 


0.162 


0.133 


0.124 


0.118 


0.107 


h 2 


20. 


0.000 


0.155 


0.143 


0.137 


0.136 


0.135 


0.133 


hz 


20. 


0.000 


0.209 


0.182 


0.171 


0.168 


0.166 


0.164 



from the data to compare it with the predictions of SUSY models, one must 
take t into account. 
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TABLE 4. The same as in the previous, but for the value of the asymmetry 
parameter A = 1.0. 











LSP 


mass 


in GeV 






Quantity 




10 


30 


50 


80 


100 


125 


250 


t 


0.0 


2.429 
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hi 


0.0 
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